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ABSTRACT 



Gulf Stream warm core rings are intense, isolated eddies which form to the north 
of the Gulf Stream from cutoff meanders which detach from the Stream. After their 
formation they are embedded in the cooler Slope Water which flows generally westward, 
impinging on the continental slope near 70° W. Observations indicate that the Slope 
Water flows westward at approximately 5 cm/s. In this study, the effect of background 
shear flows on anticyclonic eddy motion is examined using a two layer primitive 
equation model. The effect of various lateral shear profiles is considered for barotropic 
and equivalent barotropic eddies. Results indicate that eddy zonal propagation is de- 
pendent on the eddy initial vertical profile and form of the background shear lateral 
profile. Eddy meridional speeds are found to be insensitive to background shear profile 
or strength. 
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. INTRODUCTION 



The motion of western boundary current eddies, such as Gulf Stream rings, has been 
the subject of numerous observational and theoretical studies. These studies have illus- 
trated that rings move westward due to planetary rotation, have meridional motions as- 
sociated with azimuthal distortions, and are advected by background flows. 

The westward drift of an eddy due to planetary rotation or the /7-effect was discussed 
by Rossby (1948). An eddy can be viewed as a summation of Rossby waves. It decays 
due to Rossby wave dispersion. Long Rossby waves travel faster to the west than short 
waves, leading to a spreading of energy away from the eddy. Dispersion also produces 
an azimuthal distortion of the eddy which introduces a secondary propagation mech- 
anism which is northward for cyclones and southward for anticyclones. This was illus- 
trated in a series expansion by Adem (1956). This distortion was further examined in 
numerical experiments of isolated eddies by McWilliams and Flierl (1979), and Mied and 
Lindemann (1979). The effect of asymmetry on a /7-plane eddy can be seen by separat- 
ing the eddy into symmetric and asymmetric parts. Smith and Bird (1989) show this 
decomposition, which illustrates that the asymmetric part is dipolar and advects the eddy 
north or south depending on the sense of rotation. These previous studies then suggest 
that anticyclones in a quiescent background should move southwestward. Eddies em- 
bedded in a broader background flow may also have a component of propagation asso- 
ciated with mean flow advection. 

Cornillon et al. (1989) have estimated the rate of translation of anticyclonic warm 
core rings in the Slope Water north of the Gulf Stream (Figure 1 on page 2). Their es- 
timates are based on satellite images and provide ring translation rate relative to the 
surrounding Slope Water. They focus on rings which are to the east of 70° W. Rings 
and the Slope Water east of 70° W are not topographically steered by the continental 
slope as they flow to the west. Cornillon et al. (1989) also restrict their analysis to rings 
which are not interacting with the Gulf Stream. Their results, based on 1 1 rings, indicate 
a northward component (Figure 2 on page 3) of ring propagation relative to the Slope 
Water. This northward component is at odds with theoretical considerations which 
suggest that anticyclones should have a southward component of propagation. 

In another study of warm core ring translation, Auer (1987) documents the motion 
of 66 warm core rings. His study is not limited to eddies away from the Gulf Stream or 



1 




Figure 1. Gulf Stream Axis Topology: White line indicates mean Gulf Stream 

axis. Black lines indicate actual axis locations. (Courtesy of P. Cornillon) 



topography. Figure 3 on page 4 shows observations of warm core ring paths in his 
study. This figure indicates a generally southwestward warm core ring path. 

The effect of background flow on eddy motion has been studied by Matsuura and 
Yamagata (1982). Their study considers eddies which are larger than the Rossby radius 
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Figure 2. Warm Core Ring Translation \ ectors (Cornillon el al., 1989): Velocity 

(cm/s) of warm core rings relative to the Slope Water based on depth 
averaged Slope Water velocity. Solid vectors for velocities west of the 
New England Seamounts. Dashed vectors for velocities east of the 
seamounts. 

of deformation. For eddies of this size, planetary dispersion is weak and can be consid- 
ered a second order effect. They cons : oackground zonal geostrophic flow fields 

which have a corresponding merdiona a surface tilt. The surface tilt represents a 
meridional background vorticity gradient and hence introduces an additional source of 
dispersion. They then show that the weak dispersion associated with the sea surface tilt 
can exactly cancel the weak dispersion associated with the planetary /?-effect. For this 
reason, eddies in their numerical simulations with a uniform background flow do not 
disperse, but propagate as symmetric features. The applicability of this result to warm 
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Figure 3. Warm Core Ring Translation (Auer, 1987): Observed warm core and 

cold core ring tracks for 5 years. The dashed line is the 5-year mean Gulf 
Stream Landward Surface Edge (GSLSE). 



core ring propagation is questionable because of the assumption of a large eddy size. 
Warm Core rings have a radius (75 km) comparable to the Rossby radius of deforma- 
tion. 

The interaction of rings with horizontal and vertical shear was studied by Nof and 
Shi (1989) using a two-layer model based on f-plane quasigeostrophic equations. To 
obtain analytic solutions they assumed that the background shear was weak compared 
to the ring's shear. They found that eccentricity in rings can result in response to hori- 
zontal shear in the lower layer but that the rings remain circular in response to hori- 
zontal shear in the upper layer. Nof and Shi's results, applied to warm core rings, 
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suggest that the major axis of the elliptical ring would be aligned with the background 
flow if the background flow had anticyclonic shear in the lower layer. If the flow in the 
lower layer was cyclonically-sh died, the major axis would be perpendicular to the flow 
direction. Their study concentrates on eddy shape and not on eddy propagation and 
hence does not appear to explain the northward component of motion for warm core 
rings in the Slope water reported by Cornillon et al. (1989). 

In a study of the Great Red Spot on Jupiter, Ingersoll and Cuong (1981) examined 
the stability of isolated vortices in a two layer lateral shear flow. In their study, an upper 
layer vortex of size greater than the internal Rossby radius of deformation was embed- 
ded in a meridional lateral shear u{y) which is the same in both layers outside of the 
vortex. They chose a u{y) profile which was barotropically stable. An important quan- 
tity in their results is the ratio of potential to kinetic energy in the vortex. For values 

of this ratio exceeding 5.0, the vortex is unstable and is rapidly sheared apart by the basic 
PE 

flow. For ~r^=r values of approximately 1 the vortices are stable. In our study and in 
Gulf Stream rings, this value is also 0(1). None of the solutions shown here exhibit 
these instabilities. 

These various studies indicate the efFect of background shear on eddy shape, eddy 
stability, eddy decay and eddy motion. The purpose of this study is to determine the 
effects of background shear on eddy motion with emphasis on Gulf Stream warm core 
ring motion in the Slope Water. 

This paper is organized as follows. The model is described in Chapter III, with ex- 
periment numerical parameters. Numerical results obtained for initially barotropic and 
initially equivalent barotropic experiments are presented in Chapter III. Analysis of 
model results is contained in Chapters IV, V, and VI, followed by discussion and con- 
clusions. 
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II. NUMERICAL MODEL 



A. TECHNIQUE AND MODEL PARAMETERS 

Experiments are performed using a two-layer primitive equation semi-implicit nu- 
merical scheme. Motion in each layer is governed by a momentum equation: 



for layer (i= 1 upper and i=2 lower) thickness h„ transports V t , and velocities v,. The 
hydrostatic and Boussinesq approximations have been made. Density is constant in each 
layer. No mixing is allowed across the interface. Thermodynamic processes are not 
considered. Subgrid scale dissipation processes are represented by a biharmonic opera- 
tor on transport. All notation is defined in Appendix A. 

The numerical scheme (Figure 4 on page 8) has been used in numerous ocean 
mesoscale circulation studies (Hurlburt and Thompson (1980); Smith and O'Brien 
(1983)) where it is more fully discussed. It has been shown to conserve energy in the 
absence of dissipation. The model has been compared with analytic solutions for linear 
eddies (as discussed in Smith and Reid, 1982). The correct representation of Rossby 
wave dispersion processes in the model was verified. 

B. BOUNDARY AND INITIAL CONDITIONS 

A rectangular (1100 x 800 km) finite difference gridded domain is used. Grid resol- 
ution (2Ax) is 20 km. The initial state consists of a Gaussian eddy: 



in gradient balance in mid-basin. L e is the e-folding width scale for eddy. The amplitude 
of the Gaussian distribution was chosen to give a maximum velocity (v max ) of approxi- 
mately 1 m/s in all experiments. Experiments are initialized with either barotropic or 
upper layer (equivalent barotropic) eddy velocity distributions. The background shear 




[ 1 ] 



and a continuity equation: 




[ 2 ] 




[ 3 ] 
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is initially barotropic in all cases. Upper layer mean thickness is chosen to be 1000 m. 
The lower layer mean thickness is 4000 m. The first internal Rossby radius of defor- 
mation (R d ) associated with this layer thickness distribution is approximately 40 km. 
No variable bottom topography is considered. 

A radiation condition was used on the downstream (right) boundary in all model 

simulations. The radiation condition (Camerlengo and O'Brien, 1980) advects flow out 
z\jc 

of the basin (at speed -^-) when flow is outward adjacent to the boundary. The north 
and south boundaries are no-slip walls where both tangential and normal flow velocities 
are set equal to zero. 

Table 1 on page 10 provides initial conditions for the experiments. All simulations 
were integrated for a duration of 45 days. Variable parameters in the study are eddy 
vertical structure and the structure of the background shear. The biharmonic friction 
coefficient, B h , is .2jc1 0 11 . This value efficiently damps grid scale noise, leaving the eddy 
scales relatively undamped. 

The background shear flow is defined in several ways. Initial experiments (Chapter 
III C,D) have a mean flow which is barotropic (w, = u 2 ) with a linear surface tilt from 
north to south: 



h(y) = Aj- 



[4] 



where L is the domain length in the meridional direction and y is the distance from the 
southern boundary. This height field distribution has a uniform westward flow with no 
lateral shear. The flow velocity can be estimated from geostrophy: 



u = 



~g dh 
f 8y 



(10m/(s 2 ))(.4m) 

(10 _4 5 _1 )(780A:m) 



= .05 mjs 



[5] 



This velocity is consistent with that observed by Comillon et al. (1989) for the warm core 
ring region. 

The second set of experiments (Chapter III E,F) involves a lateral shear in the mean 
flow: 



, , >. 1.6 y ,2 

h(y) = — {\ - — ) 



[ 6 ] 
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Background Shear (constant, linear, or quadratic) 
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Figure 4. Model Domain for Numerical Study: Linear shear is illustrated but 

constant and quadratic shear cases are also considered. 



u = 



-g dh 

f dy 



-J'W-T )( "r ) = 



<-6g 

fL 



1 •<>&>> 

fL 1 



[7] 



The velocity shear is linear in this case with velocity equal to 5 cm/s at mid basin. As 
is discussed below, this shear' profile has anticyclonic relative vorticity, but provides no 
/7-effect (-^-) for the eddy. The lateral profile of the zonal flow is shown in Figure 5 
on page 9 Tor this case and the other shear profiles. 

The final experiments (Chapter III G,U) have a parabolic height field: 



%) = 



1.6 

3 






[S] 
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Figure 5. Velocity Profiles for Different Shear Cases: The velocity profiles for the 

constant shear, linear shear, and quadratic shear case. 



~g dh g y \2/ — 1 \ rcri 

u= — —"7 3(1 -t )( —) [9] 

For a parabolic surface distribution, the velocity distribution is quadratic. This lat- 
eral profile also has anticyclonic shear flow. In each of the sets of experiments the flow 
is initially barotropic ( u x = u 2 ) and equal to approximately 5 cm/s at mid-domain. This 
choice of velocity at mid-domain is based on the observations of Cornillon et al. (1989). 

1. Preliminary' considerations 

The characteristics of the lateral shear profiles chosen above can be examined 
to determine possible /Reflects on the eddy. In the first experiments a uniform back- 
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Table I. EXPERIMENT INITIAL PARAMETERS: Experiment names are de- 

fined by a four character code. The first two indicates barotropic (BT) or 
baroclinic (BC), the third character indicates the flow pattern (No shear, 
Constant, Linear, or Quadratic shear) and the last character indicates mid 
basin flow (W for .05 m/s and S for .108 m/s). 



Experiments 


Mean Flow Veloc- 
ity at Mid Basin 
(m/s) 


Initial Eddy Verti- 
cal Structure 


Initial Eddy Ve- 
locity (m/s) 


No Mean Flow 


BTNS 


0 


Barotropic 


1/1 


BCNS 


0 


Baroclinic 


1/0 


Constant Flow 


BTCW 


.05 


Barotropic 


1/1 


BTCS 


.108 


Barotropic 


LI 


BCCW 


.05 


Baroclinic 


1/0 


BCCS 


.108 


Baroclinic 


10 


Linear Flow 


BTLW 


.05 


Barotropic 


1/1 


BTLS 


.108 


Barotropic 


1/1 


BCLW 


.05 


Baroclinic 


1/0 


BCLS 


.108 


Baroclinic 


1/0 


Quadratic Flow 


BTQW 


.05 


Barotropic 


1/1 


BTQS 


.108 


Barotropic 


1/1 


BCQW 


.05 


Baroclinic 


1/0 


BCQS 


.108 


Baroclinic 


l/o 



ground flow is present with no lateral shear. In the second set of experiments the surface 
tilt is quadratic. In the third set of experiments the surface is parabolic. Relative 
vorticity is defined: 



dv 



dx 



du 

dy 



du 



[ 10 ] 



as no mean flow in the meridional direction is defined. Relative vorticity is equal to zero 
for the first set (constant flow) of experiments. For the second set (linear shear) relative 
vorticity is: 



10 



c = - 



[ 11 ] 



2g 

A 2 

The north-south variation of relative vorticity or a /?-efTect for the eddy, provided by the 
background flow, is zero in this case. 

In the third set experiments, the surface tilt is parabolic. The relative vorticity 
in this case is not a constant. 



du_ = 2^y_ 
dy fL 2 fL 3 



The associated vorticity gradient is: 



dj, _ 3.2 g 

8 y fL 3 



— 7.1jrl0 _13 m _1 s _1 



[ 12 ] 



[13] 



This /?-effect is associated with the background shear in the third set of experiments. 
This may augment or offset the planetary /?-effect. If this were to occur, the dispersive 
decay of the eddy due to planetary /? would be altered. This would change the eddy de- 
cay rate and direction of propagation. This was found to be the case in the results of 
Matsuura and Yamagata (1982) where the planetary dispersive decay of an eddy was 
exactly cancelled by the /TefFect associated with the shear flow. In that case the 
meridional motion of the vortex did not occur and the motion was largely zonal. In our 
case however, the shear-related /^-effect is one order of magnitude smaller than the 
planetary /^-effect. The nondispersive behavior of eddies as seen by Matsuura and 
Yamagata (1982) is not anticipated in these experiments. 
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III. MODEL RESULTS 



A. BAROTROPIC EDDY, NO SHEAR 

For comparison with subsequent experiments with background shear, a barotropic 
experiment with no mean flow (BTNS) is shown in Figure 6 on page 13. This exper- 
iment illustrates the behavior of an anticyclone as influenced by planetary and nonlinear 
self advection. As in previous isolated eddy studies, the eddy decays by Rossby wave 
radiation which spins up adjacent eddies. Long Rossby waves travel faster to west than 
do the short waves causing a distortion in the eddy. The net effect of the dispersion is 
to broaden the frontal zone in the leading edge of the eddy and to steepen the trailing 
edge. In barotropic simulations, the layers evolve in phase and the flow remains largely 
barotropic. This can be seen by comparing upper and lower layer relative vorticity, Fig- 
ure 6 on page 13. This distribution can be described by an azimuthal mode one cor- 
rection to a circular shape. While westward motion initially occurs due to /?, the 
azimuthal mode distribution leads to a southward component of propagation. By de- 
composing the eddy into an axisymmetric and nonsymmetric part and removing the 
axisymmetric part (Figure 7 on page 14), it can be seen that the nonsymmetric part as- 
sociated with azimuthal mode one is dipolar. The dipolar part then advects the sym- 
metric part southward. In this experiment the mean eddy speed is .17 km/day to the 
west and 4.57 km/day to the south. The eddy speeds for this and subsequent exper- 
iments are shown in Appendix B. The trajectory for the eddy center for this case and 
subsequent barotropic cases are shown in Figure 8 on page 15. 

Previous experiments with barotropic eddies in two layer models (McWilliams and 
Flierl (1979); Mied and Lindemann (1979)) have shown that barotropic eddies decay 
barotropically with very little energy transfer into the baroclinic mode. Barotropic 
eddies also disperse faster than baroclinic eddies. This provides greater azimuthal mode 
one distortion and hence greater southward propagation. 

B. BAROCLINIC EDDY, NO SHEAR 

A baroclinic eddy experiment with no mean flow (BCNS, Figure 9 on page 16) il- 
lustrates these differences between barotropic and baroclinic eddies. This figure shows 
that the degree of dispersion is less. That baroclinic eddies are less dispersive can be seen 
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Day 0 Upper Layer Relative Vorticity 


Day 0 Lower Layer Relative Vorticity 
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Day 22 Upper Layer Relative Vorticity 


Day 22 Lower Layer Relative Vorticity 


J 
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Day 37 Upper Layer Relative Vorticity 


Day 37 Lower Layer Relative Vorticity 



Figure 6. Barotropic Eddy, No Shear: The panels are contours of relative 

vorticity at days 0, 22, and 37. Positive values are solid, negative values 
are dashed, and values are lO^s -1 . The contour intervals equal .5x 
lCMs' 1 (day 0), .^xlO -5 ^ 1 (day 22), and .lxKH.r 1 (day 37). North is 
at the top of the Figure. 
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Figure 7. Eddy in Axisymmetric and Nonsymmetric Components (day 
22): Barotropic eddy no shear, day 22 with axisymmetric component 

removed illustrating azimuthal mode one structure. The contour is sur- 
face height anomaly after the height associated with the axisymmetric 
signature of the eddy has been removed. The contour interval is .025x 
10 _1 m. 



by considering the characteristics of Rossby waves. The Rossby wave dispersion re- 
lation is written: 



CD = 



-fik 






[14] 



where k and 1 are wavenumbers in x, y direction and R d is Rossby radius of deformation. 
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Figure 8. Trajectories for the Barotropic Cases: The symbols indicate 5 day in- 

tervals. 



For barotropic eddies the term involving the Rossby radius involves the external Rossby 
radius and is negligible compared to k and 1. This makes all the barotropic waves 
dispersive. For baroclinic eddies the internal Rossby radius is more comparable to the 
eddy scale and is not negligible. This makes a portion of the waves nondispersive. As- 
sociated with the less dispersive nature of the baroclinic eddy is less azimuthal mode one 
distortion and a smaller southward component. Westward and southward eddy speeds 
are 1.5 and 3.4 km/day respectively. The trajectory for the eddy center for this case and 
subsequent baroclinic cases are shown in Figure 10 on page 18 and Figure 11 on page 
19. Lower layer relative vorticity (Figure 9 on page 16) illustrates that lower layer cir- 
culations are the result of vortex stretching of quiescent lower layer fluid by the moving 
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Day 0 Upper Layer Relative Vorticity Day 0 Lower Layer Relative Vorticity 





Day 22 Upper Layer Relative Vorticity Day 22 Lower Layer Relative Vorticity 





Day 37 Upper Layer Relative Vorticity Day 37 Lower Layer Relative Vorticity 

Figure 9. Baroclinic Eddy, No Shear: The panels are contours of relative 
vorticity at days 0, 22, and 37. Positive values are solid, negative values 
are dashed, and values are 10~ s s _1 . The contour intervals equal .5x 
ICKr 1 (day 0), .125xI0- 5 5 _1 (day 22), and .lxl0~ 5 ,s _1 (day 37). North is 
at the top of the figure. 
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upper layer eddy. As argued by Cushman- Roisin et al. (1989), fluid columns to the west 
of the eddy are vortex squashed as the anticyclone goes over them. Fluid columns ini- 
tia under the anticyclone are vorte stretched as the anticyclone moves westward. 
Custiman-Roisin et al. (1989) argue that these secondary circulations combine to advect 
the upper layer anticyclone to the south (see Figure 12 on page 20). 

A comparison with the reduced gravity results of Smith and Reid (1982) where upper 
layer eddies over a motionless lower layer are considered indicates that this lower layer 
effect does contribute to the southward motion. In Smith and Reid (1982) a comparable 
size and strength eddy moved in the meridional direction at 1 km/ day in contrast to 3.4 
km/day in BCNS. 

C. BAROTROPIC EDDY, CONSTANT SHEAR 

For comparison with the lateral shear flow cases below, several experiments with 
constant .05 m/s background flow are discussed. Figure 13 on page 21 shows relative 
vorticity for barotropic case BTCW. As in case BTNS, rapid barotropic dispersion gives 
rise to a trailing cyclonic vortex to the northeast of the anticyclone. The anticyclone 
moves south (4.89 km/day) and west (4.02 km/day). 

The meridional propagation speed of the eddy is comparable to no mean flow case 
BTNS. The zonal speed is substantially higher than in BTNS. The eddy speed to the 
west (4.02 km/day) is 10% less than the mean flow speed (5 cm/s) (5 cm/s equals 4.32 
km/day) plus zonal eddy speed in the no shear case BTNS (.17 km/day) suggesting that 
the two efFects on eddy motion are not simply additive. Appendix B shows the zonal 
and meridional speeds For this and other barotropic cases. 

The surFace height anomaly Fields For this experiment show that a divergent region 
appears to the west oF the eddy. The eddy thus has a substantial efFect on the back- 
ground flow upstream and downstream oF the eddy. As in experiment BTNS, the layers 
evolve in phase with no significant transfer of energy into the baroclinic component. 

D. BAROCLINIC EDDY, CONSTANT SHEAR 

Figure 14 on page 22 shows the evolution of corresponding baroclinic case BCCW 
For constant background flow (.05 m/s). As with no mean flow (baroclinic case BCNS), 
the eddy is less dispersive than the preceding barotropic case. A weak trailing cyclone 
appears to the east of the upper layer anticyclone in relative vorticity. This cyclone is 
also evident as a northward distortion in the mean flow lines to the east of the 
anticyclone in the surFace height anomaly. This plot also shows that while an upstream 
efFect on the mean flow is caused by the eddy, no downstream efFect (as in BTCS) is 
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Figure 10. Trajectories for Baroclinic Cases 1-4: The symbols indicate 5 day in- 

tervals. 



seen. Like the barotropic cases, meridional speed is unaltered by the presence of the 
shear (3.27 km/day), but the zonal speed (6.21 km/day) is a factor of four greater than 
in the no background flow case BCNS. As in the barotropic constant shear case 
(BTCW), the zonal speed in this case is less (5% in this case) than the sum of the pre- 
vious westward component in BCNS (1.5 km/day) and the constant 5 cm/s (4.32 
km/day) background flow'. 

The dipolar lower layer flow' previously seen in experiment BCNS (Figure 9 on page 
16) is very comparable in this experiment. Figure 15 on page 23 show's velocity vectors 
in the lower layer for this case. The cyclonic flow' under the eddy is approximately 8 
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Figure 11. Trajectories for Baroclinic Cases 5-7: The symbols indicate 5 day in- 

tervals. 



cm/s. The anticyclonic flow is weaker (approximately 2 cm/s). Despite the southward 
advective effect associated with this lower dipole, this baroclinic anticyclone moves to 
the south slower than the previous barotropic case (BTCW). The faster dispersion in the 
barotropic mode explains this difference. 

E. BAROTROPIC EDDY, LINEAR SHEAR 

Barotropic eddy motion in a linearly sheared background flow (Figure 16 on page 
24) is very comparable to that in the constant shear case BTCW. This can be seen by 
comparing trajectories in Figure 8 on page 15. Meridional speeds arc nearly the same 
with zonal speed slightly higher at 5.5 km/day (vs. 4.0 km/day in the constant shear 
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Figure 12. Vortex Stretching of Quiescent Loner Layer Fluids (Cushman- Roisin 
et al., 1989): Impact of the westward migration of an upper-layer eddy 

on the lower layer and the reaction of the induced lower-layer relative 
vorticity on the eddy drift: squeezing and stretching under (a) an 
anticyclone and (b) a cyclone. 



case). The higher zonal velocity is associated with the increase in background flow ve- 
locities as the eddy moves southward. The divergent region upstream previously seen in 
the constant shear case is also seen in this experiment as is the trailing cyclonic vortex. 
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Day 0 Upper Layer Relative Vorticity 




Day 22 Upper Layer Relative Vorticity 




Day 37 Upper Layer Relative Vorticity 





Day 22 Surface Height 




Figure 13. Barotropic Eddy, Constant Shear: The left panels are contours of rel- 

ative vorticity. Positive values are solid, negative values are dashed, 
and values are The contour intervals are .5x1 O^s -1 (day 0), .25x 

10 ~ s s~ l (day 22), and .lxlO -5 ^ 1 (day 37). The right panels are contours 
of surface height (cm). The contour intervals are 5 cm (day 0), 2.5 cm 
(day 22), and 2.5 cm (day 37). North is at the top of the figure. 
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Day 37 Upper Layer Relative Vorticity Day 37 Surface Height 

Figure 14. Baroclinic Eddy, Constant Shear: The left panels are contours of rel- 

ative vorticity. Positive values are solid, negative values are dashed, 
and values are KHs -1 . The contour intervals are .5xl0 _5 5 _1 (day 0), .25x 
10- 5 5 _1 (day 22), and .lxKTV 1 (day 37). The right panels are contours 
of surface height (cm). The contour intervals are 5 cm (day 0), 5 cm 
(day 22), and 5 cm (day 37). North is at the top of the figure. 
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Figure 15. Lower Layer Velocity Vectors and Isotachs for Baroclinic No Shear (day 
22): Contours interval is 1 cm/s with arrows indicating magnitude and 

direction. 

F. BAROCLINIC EDDY, LINEAR SHEAR 

As in the preceding section, a linear shear (Figure 17 on page 25) only slightly 
changes meridional baroclinic eddy propagation from the constant (no shear) case. The 
trajectory (Figure 10 on page 18 and Figure 1 1 on page 19) for this case is comparable 
to the no shear case, with zonal speed slightly higher (6.9 vs. 6.2 km/day). These speeds 
to the west are higher than the corresponding barotropic cases. This is true in the 
quadratic shear cases also. The reasons for the greater speeds in baroclinic cases are 
discussed below. 
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Day 37 Upper Layer Relative Vorticity 



Day 37 Surface Height 



Figure 16. Barotropic Eddy, Linear Shear: The left panels are contours of rela- 

tive vorticity. Positive values are solid, negative values are dashed, and 
values are lO^s -1 . The contour intervals are JxlO 5 ^ 1 (day 0), .125x 
lO -5 ^ -1 (day 22), and ^xlO -5 ^ 1 (day 37). The right panels are contours 
of surface height (cm). The contour intervals are 5 cm (day 0), 2.5 cm 
(day 22), and 2.5 cm (day 37). North is at the top of the figure. 
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Day 37 Upper Layer Relative Vorticity 


Day 37 Surface Height 



Figure 17. Baroclinic Eddy, Linear Shear: The left panels are contours of relative 

vorticity. Positive values are solid, negative values are dashed, and 
values are lO -5 ^ -1 . The contour intervals are .5x1 0'V 1 day 0), .25x 
KUV 1 (day 22), and . Ixl0 -5 5 -1 (day 37). The right panels are contours 
of surface height (cm). The contour intervals are 5 cm (day 0), 5 cm 
(day 22), and 5 cm (day 37). North is at the top of the figure. 
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G. BAROTROPIC EDDY, QUADRATIC SHEAR 

The presence of quadratic shear leads to a slightly higher meridional eddy speed over 
previous barotropic cases. Zonal speed averages 7.8 km/day (40% greater than the cor- 
responding linear shear case). A comparison of Figure 18 on page 27 with the linear 
shear case (Figure 16 on page 24) shows only minor differences in the spatial eddy 
structure. 

Figure 19 on page 28 shows velocity vectors and isotachs for this case. A compar- 
ison of upper and lower layers indicate that the flow remains barotropic during the ex- 
periment. Velocities in the divergent region to the west of the eddy approach zero with 
time. On day 22, the strongest flow is on the trailing side of the eddy. At this time the 
eddy center is located at the latitude of 10 cm/s zonal flow. 

H. BAROCLINIC EDDY, QUADRATIC SHEAR 

As in previous sections, meridional speed (3.2 km/day) is unaltered by the form of 
the shear and a comparison of Figure 20 on page 29 with the linear shear case 
(Figure 17 on page 25) shows only minor differences. Zonal speed (8 km/day) is higher 
than the linear shear case. This is consistent with the stronger flows near the southern 
boundary 7 (Figure 21 on page 30) in this case relative to the previous shear flows. In 
contrast to previous cases this speed is also comparable to the barotropic quadratic 
shear case (BTQW). In previous cases the zonal speed was higher in baroclinic cases. 

To examine how background flows modify the eddy motion, the trajectories for 
baroclinic and barotropic cases were recomputed with the initial background shear flow 
removed. Figure 22 on page 31 and Figure 23 on page 32 show these trajectories for the 
barotropic and baroclinic cases respectively. These figures show that a constant shear 
nearly adds linearly to eddy propagation. The trajectories for linear and quadratic shear 
cases indicate an eastward component of motion induced in the eddies by background 
flow. For quadratic shear barotropic case BTQW, for example, the eddy is displaced 120 
km in 45 days or an eastward speed of 2.7 km/day relative to the background westward 
flow. It is interesting to compare the trajectory for the quadratic shear baroclinic case 
BCQW and the corresponding barotropic case BTQW. While the meridional component 
of motion in BTQW is higher the zonal components are nearly identical. This is also 
supported by the mean zonal components for these cases (7.8 amd 8.0 km/day) in 
Table 2 on page 48 and Table 3 on page 48. 
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Day 0 Upper Layer Relative Vorticity 




Day 22 Upper Layer Relative Vorticity 




Day 0 Surface Height 





Day 37 Upper Layer Relative Vorticity 



Day 37 Surface Height 



Figure 18. Barotropic Eddy, Quadratic Shear: The left panels are contours of 

relative vorticity. Positive values are solid, negative values are dashed, 
and values are 10- 5 s -1 . The contour intervals are .5xl0‘ 5 s -1 (day 0), .25x 
10- 5 5-* (day 22), and .IxlCKr 1 (day 37). The right panels are contours 
of surface height (cm). The contour intervals are 5 cm (day 0), 5 cm 
(day 22), and 5 cm (day 37). North is at the top of the figure. 
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Shear (day 22): Contour interval is 5 cm/s with arrows indicating 

magnitude and direction. 
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Day 22 Upper Layer Relative Vorticity 




Day 37 Upper Layer Relative Vorticity 



Day 22 Surface Height 




Figure 20. Baroclinic Eddy, Quadratic Shear: The left panels are contours of 

relative vorticity. Positive values are solid, negative values are dashed, 
and values are lO -5 ^ 1 . The contour intervals are .5xl0 _5 5 _1 (day 0), .25x 
KKs- 1 (day 22), and .1x10- 5 5 _1 (day 37). The right panels are contours 
of surface height (cm). The contour intervals are 5 cm (day 0), 5 cm 
(day 22), and 5 cm (day 37). North is at the top of the figure. 
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Figure 21. Lower Layer Velocity Vectors and Isotachs for Baroclinic Quadratic 
Shear (day 22): Contours interval is 1 cm/s with arrows indicating 

magnitude and direction. 
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Figure 22. Trajectories for Barotropic Cases Minus Background Shear: The tra- 

jectories for barotropic cases were recomputed with the initial back- 
ground shear flow removed. The symbols indicate 5 day intervals. 



31 




Figure 23. Trajectories for Baroclinic Cases Minus Background Shear: The tra- 

jectories for baroclinic cases were recomputed with the initial back- 
ground shear flow removed. The symbols indicate 5 day intervals. 
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IV. AZIMUTHAL MODE ANALYSIS 



The zonal and meridional average eddy speeds (Table 2 on page 48 and Table 3 on 
page 48) are consistent with azimuthal distortions in the eddies. In Figure 24 on page 
35 and Figure 25 on page 36 the eddies have been decomposed into axisymmetric and 
asymmetric components. Only the asymmetric portion of the eddy is shown in these 
Figures. The dipolar structure in the majority of the experiments indicates that an 
azimuthal mode one structure is present. As discussed above, this distortion is caused 
by Rossby wave dispersion. The /?-e fleet associated with quadratic shear has little effect 
in offsetting this dispersive behavior. The orientation of the dipolar part provides 
southward and westward eddy motion in most baroclinic cases. A comparison of the 
orientation of the dipolar structures in Figure 25 on page 36 indicates that the back- 
ground shear has little effect in rotating the azimuthal mode one structure. The higher 
westward eddy speeds for the linear and quadratic shear cases (Table 3 on page 48) thus 
are related to the higher background flow velocities and not to a rotation of the mode 
one structure by the shear. 

Azimuthal mode one structure is also present in the barotropic cases. Experiments 
BTNS and BTCW (Figure 24 on page 35) have dipolar corrections to a circular sym- 
metric shape. The orientation of the dipole gives southeastward motion in the case of 
BTNS and southward for BTCW. The trajectory for BTNS (Figure 8 on page 15) shows 
eastward motion near the end of the run. In contrast to baroclinic cases (all show 
azimuthal mode one structure), barotropic cases BTLW and BTQW have higher 
azimuthal mode structure. The presence of two highs and two lows in Figure 24 on page 
35 suggests that these eddies had a more elliptical distribution. The elliptical distrib- 
utions do not contribute to eddy motion as does azimuthal mode one distribution. A 
comparison of westward and southward mean speeds (Table 2 on page 48 and Table 3 
on page 48) for barotropic and baroclinic eddies reveals that barotropic eddies move 
faster to the south and slower to the west than baroclinic eddies. This seems inconsistent 
with the fact the barotropic eddies which propagate further south than baroclinic eddies 
are advected westward by higher background velocities. The enhanced westward speed 
for baroclinic eddies appears related to the azimuthal mode orientation which augments 
the westward /^-induced speed. The orientation of the azimuthal mode one structure in 
the barotropic cases does not give a westward component of motion. In addition, it was 
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demonstrated above that barotropic eddies in linear or quadratic shear have eastward 
components of motion relative to the background flow. 
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BTQW day 22 

Figure 24. Azimuthal Mode Decomposition for Barotropic Cases (day 22): The 

contour is surface height anomaly after the height associated with the 
axisymmetric signature of the eddy has been removed. The contour 
interval equal .025xl0 _1 m for the Barotropic No Shear (BCNS) case, 
.020xl0 _1 m for the Barotropic Constant Weak (BTCW) case, 
.025xl0 _1 m for the Barotropic Linear Weak (BTLW) case, and .025x 
ICHm for the Barotropic Quadratic Weak (BTQW) case. Weak indi- 
cates mid basin flow of .5 m/s. North is at the top of the figure. 
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BCLW day 22 BCQW day 22 

Figure 25. Azimuthal Mode Decomposition for Baroclinic Cases (day 22): The 

contour is interface height anomaly after the height associated with the 
axisymmetric signature of the interface has been removed. The contour 
interval equal 2.5xl0~ 1 m for the Baroclinic No Shear (BCNS) case, 2.5x 
lO-’m for the Baroclinic Constant Weak (BCCW) case, 2.5xl0-'m for 
the Baroclinic Linear Weak (BCLW) case, and 2.5xlO _1 m for the 
Baroclinic Quadratic Weak (BCQW) case. Weak indicates mid basin 
flow of .5 m/s. North is at the top of the figure. 
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V. TIME DEPENDENT EDDY MOTION 



A. BAROTROPIC CASES 

Figure 26 on page 38 shows the zonal eddy speeds as a function of time for all 
barotropic cases. This figure shows that the eddy speed is time dependent. The eddy 
speeds plotted in this figure have been averaged into 5 day intervals. In the absence of 
any background flow (case BTNS), the highest eddy speeds are during the first ten days. 
An increase in speed during the first 10 days is seen in all cases. This is associated with 
an adjustment period during which the eddy adjusts from the initial axisymmetric shape 
to a dispersing eddy. Figure 26 on page 38 also shows that the effect of a constant 
background flow (case BTCW) is to almost add linearly (within 10%) to eddy speed in 
the absence of shear. The presence of linear and quadratic shear causes eddy speed to 
increase with time after day 25 as the eddy propagates into the higher velocity region. 

The corresponding meridional motion for all barotropic cases is shown in 
Figure 27 on page 39. Meridional speeds are also time dependent but little difference 
exists between the various cases. Eddy adjustment during the first 15 days causes 
meridional speed to increase during this period. Meridional speed decreases after day 15 
associated with rapid barotropic dispersion. 

These zonal and meridional speed changes are the result of several factors. Initially, 
eddy shape adjustments lead to increases in both zonal and meridional speed. Subse- 
quently dispersive weakening of the eddy over time decreases both zonal and meridional 
speeds. Finally an increase in westward speed is associated with increased zonal back- 
ground flow (in the linear and quadratic shear cases BTLW and BTQW). For the zonal 
velocity then, initially dispersive adjustment and decay is important and later, stronger 
background flow magnitude is important. 

B. BAROCLINIC CASES 

The time dependent eddy speeds for all baroclinic cases are shown in Figure 28 on 
page 40 (zonal) and Figure 29 on page 41 (meridional). The baroclinic eddy speeds also 
show large changes in speed with time. The baroclinic zonal speeds also show a large 
initial adjustment period during the first 15 days. In contrast to the barotropic cases 
where eddy speeds increase during the adjustment period, baroclinic cases BCCW and 
BCLW show decrease in zonal eddy speed initially. Beyond day 15 zonal speeds are rel- 
atively constant. The less dispersive baroclinic eddies do not show the slowing seen in 
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Figure 26. Zonal Time Dependent Speeds for Barotropic Eddies: Mid basin flow 

is .05 m/s. Positive values represent North for meridional flows and 
East for zonal flows. 



the barotropic cases. There is a slight tendency for baroclinic zonal speeds to increase 
near the end of the of the 40 day period. All baroclinic meridional speeds decrease to less 
than several km/day during the first 10 days. This is in contrast to barotropic meridional 
speeds which increased during the first 20 days. Maximum baroclinic southward speeds 
reach 4 km'day in contrast to 8 km/day for barotropic cases. 

C. TIME AVERAGED ZONAL AND MERIDIONAL EDDY SPEED 

Figure 30 on page 42 summarizes average meridional and zonal eddy speeds for 
barotropic and Figure 31 on page 43 for baroclinic cases. These figures illustrate that 
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Figure 27. Meridional Time Dependent Speeds for Barotropic Eddies: For .05 m/s 

mid basin flow 



eddy meridional speed is relatively insensitive to the form of background shear, but zonal 
speed is dependent on the shear profile. 

D. COMPARISON WITH OTHER MODEL RESULTS 

The model of Cushman-Roisin et al. (1989) suggests several eddy propagation 
tendencies with which these results can be compared. Cushman-Roisin et al. (1989) ar- 
gues that eddy propagation can be influenced by two mechanisms which depend on the 
eddy vertical structure. For an anticyclone in the upper ocean only, vortex squashing 
of fluid columns in the lower layer will create anticyclonic vorticity under the leading 
edge of the eddy. Fluid columns in the lee of the eddy will acquire cyclonic vorticity 
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Figure 28. Zonal Time Dependent Speeds for Baroclinic Eddies: For .05 m/s mid 

basin flow 



through vortex stretching. These tendencies were illustrated in Figure 12 on page 20 
(from Cushman- Roisin et al. 1989). The combined effect is to advect the eddy to the 
south. For barotropic eddies the fluid columns must move north and south around the 
eddy in each layer. Fluid columns moving northward acquire anticyclonic vorticity. 
Fluids columns moving southward around the eddy acquire cyclonic vorticity as the 
Coriolis parameter decreases. The combined effect of these tendencies is to augment the 
westward speed of the barotropic anticyclones. 

For no background flow, the results of baroclinic experiment BCNS illustrates the 
Cushman-Roisin et al. (1989) baroclinic argument (see Figure 9 on page 16). In the 
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Figure 29. Meridional Time Dependent Speeds for Baroclinic Eddies: For .05 m s 

mid basin flow 



presence of background flows, these tendencies still apply. Fluid columns in the lower 
layer should acquire the same vorticity tendencies as they go under the eddy, as the up- 
per layer eddy moves relative to the lower layer mean fluid speed. The Cushman- Roisin 
et al (1989) argument is thus unmodified by a barotropic background flow. The 
meridional speed for baroclinic cases is nearly independent of background shear flow 
magnitude and lateral profile. 

The Cushman-Roisin et al. (1989) barotropic argument can also be modified for 
background flow. For barotropic background flow with speeds greater than the eddy 
speed, fluid columns overtaking the eddy must go around the eddy to the north and 
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Figure 30. 



Summary of Time Averaged Meridional and Zonal Speeds for Barotropic 
Eddies. 



south. The northward (southward) displaced columns acquire anticyclonic (cyclonic) 
vorticity. These effects augment the westward ^-induced eddy speed. Figure 6 on page 
13 (BTNS) does not show any indication of these vorticity tendencies. Furthermore it 
was shown in Figure 22 on page 31 that barotropic anticyclones had an eastward com- 
ponent of motion relative to the background flow. 

E. COMPARISON WITH OBSERVATIONS 

For comparison with the observations of warm core ring propagation by Cornillon 
et al. (1989) and Auer (1987) average speed and direction for barotropic and baroclinic 
eddy experiments were computed. The averages were based on the experiments with no 
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Experiments 



Figure 31. Summary of Time Averaged Meridional and Zonal Speeds for Baroclinic 
Eddies. 

shear or shears with 5 cm/s at mid basin. The resulting speeds are 6.7 km/day at 230° 
from north for barotropic eddies and 6.4 km/day at 241° for baroclinic eddies. Based on 
observations of 66 warm core rings, Auer (1 S7) found a mean speed of 2.4 km/day at 
256°. 

Eddy speeds in experiments here are considerably higher than the observations sug- 
gest. Higher speeds in these experiments suggest that background in the model flow ex- 
ceed those to the north of the Gulf Stream. 

Results here do not explain the northward component of motion for warm core rings 
east of 70° W observed by Comillon el al . (1989). Their study, based on 11 warm core 
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rings, gives a mean 8.5 km/day at 282° from north illustrating northward motion of rings 
relative to the slope water. None of the shear flows chosen here provide this tendency. 
Observations of warm core ring velocities from other investigators are shown in Ap- 
pendix C. This table illustrates that estimates of warm core ring speeds vary widely 
ranging 5 to 7 km/day. 
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VI. SUMMARY 



The effect of background shear flows on ocean eddy motion has been considered 
using a two-layer numerical model. The results indicate that eddy meridional motion is 
not dependent on the magnitude or form of the background shear profile. This is true 
for barotropic and baroclinic experiments. Initially barotropic eddies have higher 
meridional speeds than baroclinic eddies due to faster Rossby wave dispersion in 
barotropic eddies. 

Eddy zonal speed depends on the initial vertical velocity profile of the eddy. 
Barotropic eddies propagate slower to the west in all cases than do baroclinic eddies 
despite the fact that barotropic eddies propagate further into the high velocity regions 
of the background flow. The higher baroclinic speeds are associated with azimuthal 
mode one distortions. These distortions have an orientation which advects the baroclinic 
eddy southwestward. The /? induced westward motion is thus augmented by this 
azimuthal mode one distortion. Barotropic eddies can have azimuthal mode one dis- 
tortions but the orientation of these does not give westward motion in these exper- 
iments. 

None of the results in this thesis appear to explain the observations by Cornillon et 
al. (1989) of warm core rings propagating northwestward relative to the background 
Slope Water. The eddies in that study were presumably not interacting with the Gulf 
Stream or continental slope topography. The close proximity of the Gulf Stream with 
the continental slope however makes these interactions likely. From this standpoint, the 
results of Stern and Flierl (1987) may be applicable. Their results suggest northward 
motion for anticyclones interacting with the cyclonic shear on the north side of the Gulf 
Stream. They showed that these interactions could induce northward eddy motion for 
eddies as far away from the Gulf Stream as several Rossby radii (100-150 km). 
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APPENDIX A. SYMBOLS AND NOTATION 



B, 


biharmonic lateral friction coefficient = .2xl0‘ 1 


fo 


Coriolis parameter for mean latitude 40°N 


g 


Gravitational acceleration 




g' 


Reduced gravitational acceleration 


- g P2 /' - -02 


K 


Instantaneous upper (i= 1) and lower (i= 2) layer thickness 


H, 


Upper (i= 1) and lower (i= 2) layer 


mean thickness 


L e 


e-folding scale for the eddy 


= 40 km 


L 


North-south domain length 


= 780 km 


Pi 


Pressure in the upper layer 


= g(/2, + h 2 + d) 


Pi 


Pressure in the lower layer 


-sf 

No 

1 

II 


K 


First internal Rossby radius of deformation 






-Jh 1 + h 2 =40km 


K 


Rossby number 


ii 


v , 


Velocities in the x and y directions 




u„ v t 


Transport in the x and y directions 
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x,y 



y 

Ax, Ay 
A t 



P, 



V 

V 4 



Cartesian coordinates directed N and W respectively 
nondimensional eddy size 



Grid spatial resolution 
Model time step 
Density in i th layer 
Gradient operator 
biharmonic operator 



(2Ax = 20 km) 
4200 sec 



dx dy 

a 4 e 4 

dx 4 dy 4 



47 



APPENDIX B. TIME-AVERAGED EDDY TRANSLATIONAL SPEEDS 



Table 2. INITIALLY BAROTROPIC ANTICYCLONES 



Exp 


Zonal 

Speed 

(km/day) 


Meridional 

Speed 

(km/day) 


Average 

Speed 

(knYday) 


Lateral Shear 
(m/sec) 


BTNS 


-0.17 


-4.57 


5.53 


No Shear 


BTCW 


-4.02 


-4.89 


7.80 


Constant (.054) 


BTCS 


-8.89 


-5.31 


12.12 


Constant (.108) 


BTLW 


-5.50 


-5.12 


9.11 


Linear (.054) 


BTLS 


-11.12 


-5.66 


14.67 


Linear (.108) 


BTQW 


-7.80 


-5.77 


10.88 


Quadratic (.054) 


BTQS 


-14.83 


-5.43 


18.15 


Quadratic (.108) 



Table 3. INITIALLY BAROCLINIC ANTICYCLONES 



Exp 


Zonal 

Speed 

(km/day) 


Meridional 

Speed 

(km/day) 


Average 

Speed 

(km/day) 


Lateral Shear 
(m/sec) 


BCXS 


-1.51 


-3.38 


7.17 


No Shear 


BCCW 


-6.21 


-3.27 


12.54 


Constant (.054) 


BCCS 


-11.26 


-3.33 


17.97 


Constant (.108) 


BCLW 


-6.94 


-3.19 


13.32 


Linear (.054) 


BCLS 


-12.55 


-3.13 


18.84 


Linear (.108) 


BCQW 


-8.00 


-3.15 


13.93 


Quadratic (.054) 


BCQS 


-14.69 


-3.18 


19.18 


Quadratic (.108) 
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APPENDIX C. WARM CORE RING SPEEDS 



Table 4. OBSERVED WARM CORE RING SPEEDS 



Investigator 


Observed Speed 
(km/day) 


Mean Speed 
(km/day) 


Number of Rings 


Cheney (1976) 


- 


5.0 


1 


Bisagni (1976) 


6.1 - 8.3 


7.0 


13 


Lai and 

Richardson (1977) 


3 - 7 


- 


20 


Halliwcll and 
Mooers (1979) 


1.7 - 13.0 


5.2 


14 


Brown et al. (1886) 


1.4 - 11.9 


5.6 


87 


Auer (1987) 


- 


2.4 


66 


Cornillon et al. 
(1989) 


- 


8.5 


11 
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